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Following the subtraction procedure for manifolds with boundaries, we calculate by variational 
methods, the Schwarzschild-de Sitter and the de Sitter space energy difference. By computing the 
one loop approximation for TT tensors we discover the existence of an unstable mode even for the 
non-degenerate case. This result seems to be in agreement with the sub-maximal black hole pair 
creation of Bousso-Hawking. The instability can be eliminated by the boundary reduction method. 
Implications on a foam-like space are discussed. 
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, I. INTRODUCTION 

o . 

^v^j ■ An intriguing property of quantum physics is the particle creation generated by external fields, like a constant 
electric or magnetic field, or by quantum fluctuation of the vacuum. In this last case only virtual particles are 
involved. However when the energy scale is large enough virtual particles can be transformed into real. In line of 
principle the same mechanism can be shared by the gravitational field where virtual black holes can be created 
\ and annihilated in analogy with particle physics. This particular phenomenon has been investigated in different 
contexts and in particular when a cosmological constant is introduced . In this example the process is mediated 
by the corresponding gravitational instanton, and the semiclassical nucleation rate for a pair on a given background 
C*~) ■ is given by 
> ' 

00 ! r = Aexp [- (I inst - I back )} . (1) 

hnst is the classical action of the gravitational instanton mediating the pair creation, Iback is the action of the 
background field, and A is the prefactor containing quantum corrections. For the de Sitter (dS) space the quantum 
creation of black holes leads to the discovery of an unstable mode in the physical sector, when one- loop approximation 
is considered ||-^|. This quantum instability is related to the S 2 x S 2 instanton responsible for the pair creation 
process. This instanton, termed the Nariai instanton B, is nothing but the extreme Schwarzschild-de Sitter (SdS) 
solution written in another system of coordinates. This instability leads to spontaneous nucleation of black holes 
signaling a transition from a false vacuum to a true one S . This transition is possible when the energy stored in the 
boundaries is the same for both spaces |)| . However as remarked in Ref . Q , the nucleation appears with a temperature 

Tpair = y ^ different from the temperature of the heat bath, which is the dS space with Tds = 5^y4f ■ This does 
not happen, for example, in the hot Minkowski space where the nucleated black hole has the same temperature as the 
/S 1 heat bath 1(1. The same situation holds even when we consider a negative cosmological constant, i.e. Anti-de Sitter 
(AdS) space. In fact to spontaneously nucleate a black hole, which has an intrinsic temperature Ts-AdSi the same 
temperature has to be imposed to the AdS space fTl| [l3|] • However in Ref. jl4| , we have shown that a semi-classical 
instability (WKB) appears for Minkowski space even at zero temperature, provided that boundary conditions be 
energy preserving. The same semi-classical instability appears also for the AdS space ]l5|] with the same energy 
condition. An interesting common feature between these cases comes from the expression of the mixed Ricci tensor 
Rf computed in these different backgrounds, namely the Schwarzschild, S-AdS and SdS metrics respectively. Indeed 
in the first case, Rf has components: 

(_2MGMGMG\ 

while the case with the cosmological constant is 

where the upper case and the lower case are related to the SdS and the S-AdS metrics respectively and b 2 = 3/A c . 
A c is the positive cosmological constant. It is straightforward to note that the only difference between these different 
tensors is in the presence of the cosmological term. This means that the source of instability residing in the first 



1 



component appears even in the other cases. However, only at one loop is possible to reveal the presence or the absence 
of such an instability. We recall that one loop computations of the energy in this context represent a Casimir-like 
energy which measures vacuum fluctuations. Following Refs. |L4|,[n|, we will consider a constant time slice £ of the 
SdS manifold whose perturbations at £ in absence of matter fields define quantum fluctuations of the Einstein- 
Rosen bridge. Indeed, as will follow in Section || even though the SdS is not asymptotically flat, the hypersurface 
£ defines a wormhole with topology S* 2 x /, where / C R is a sub-interval of R. This is a consequence of having a 
cosmological radius which sets an upper bound to the radial coordinate. To this purpose we will fix our attention on 
a Hamiltonian with boundary 

H T = H^ + H as = d 3 x{NH+N t H l ) + ff as , (4) 
where N is called the lapse function, Ni is the shift function and 

= ~H 

Hon represents the energy stored into the boundary. The aim of this paper is the evaluation of 



E 



- ls (M, b) = E dS (b) + AE S J S (M, b) lclassical + AE S J S (M, b)^ loop , (6) 



representing the total energy computed to one-loop in a SdS background. E ds (b) is the reference space energy, i.e. 
the de Sitter space. A-Ejg (M, b)^ classical is the energy difference between the SdS and the dS metrics, stored in 
the boundaries and AE^ S (M 7 b)^_ loop is the quantum correction to the classical term. The rest of the paper is 
structured as follows, in section [n] we compute the quasilocal energy for the SdS space, in section III we give some of 



the basic rules to perform the functional integration to evaluate the energy density of the Hamiltonian approximated 



up to second order in the SdS background, in section IV we look for stable modes of the spin-two operator acting on 
transverse traceless tensors, in section |y| we show the existence of only one negative mode under suitable conditions 



and we compute the energy density for stable modes, in section VI, we confirm the existence of one negative mode 



for the extreme SdS background, namely the Nariai metric, als o in our framework and we give a computation for the 



stable part in analogy with its non-extreme sector, in section VII we fi nd a critical radius below which we have a 



stabilization of the system. We summarize and conclude in section VIII 



II. QUASILOCAL ENERGY FOR THE SDS SPACE 

In this section we fix our attention to the classical part of Eq. @ . We begin to define the line element 

ds 2 = -f (r) dt 2 + f (r)" 1 dr 2 + r 2 dtf, (7) 



referred to the SdS metric, where 



b 2 

For A c = the metric describes the Schwarzschild metric, while for M = 0, we obtain the de Sitter metric (dS) 



/M-|i-^-£|. («) 



ds 2 = -(l--)dt 2 +(l-^ dr 2 + rW. (9) 



b 2 J V b2 J 

The gravitational potential / (r) admits three real roots. One is negative and it is located at 

r _ = -^6ooB(— J, (10) 



x In Appendix we will report the details concerning the Kruskal-Szekeres description of the SdS manifold. 
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2 , {9 + 4ir\ 2 , . . . 

r + = — ^ocos — - — I , r ++ = —=bcos — I (11) 



while 

= V3"^\ 3 V3"^V3 
are associated to the black hole and cosmological horizons respectively, with 

cos = -UIGVs/b. (12) 

However in the wormhole language, we will say that r + is the inner throat and r ++ is the outer throat. Note also 
that the hypersurface £ is described by the three-dimensional wormhole whose metric is 

ds 2 = /(r) _1 dr 2 +rW, (13) 

where / (r) is given by Eq. (|J) . A relation between the three roots is given by 

b 2 = r\ + r + r ++ + r\ + 
2Ml 2 p b 2 = {r+r ++ ){r++r ++ ) . (14) 
= r_ + r + + r ++ 



Thus , we can write 



with 



f M = ( r ~ r +) ( r ~ r + + ) ( r + 7 '+ + r ++) ' ( 15 ) 



r a- < r < r_ 



7T 37T 

2'T 



Since r + is a monotonic increasing function of 9, while r ++ is a monotonic decreasing with 

r l +e[o,b] 

r l ++ e[0,b}' 



(16) 



(17) 



in order to have the inequality ( |l^ ) preserved, we have to consider 9 £ [? , 7r] . Indeed when S [7T, ^1 the inequality 
( |l6| ) is reversed and the meaning of the internal and external roots is exchanged. Thus the de Sitter region delimited 
by the bound of Eq. (|l6|) , at time fixed, can be represented as in Fig.|l|. 




r=r + w_^ r=r + 
r=r ++ 

FIG. 1. The geometry of the constant time slice embedded in flat space with a polar angle suppressed. Isometric copies of 
this surface can be smoothly joined at the throats, producing a periodic S 2 x R 1 spatial topology. 

A common value is reached when 9 = it where r + = r ++ = b/y/3 and the metric is termed extreme. This particular 
case will be discussed in section |Vl]. However, instead of looking at the de Sitter region with the topology of Fig|l], 
we will look at the Einstein-Rosen bridge corresponding to the inner bifurcation surface depicted in Fig.0. 




r=r++ r=r 

FIG. 2. The same representation of Figjl| but with the cosmological "wormhole mouths'" placed in antipodal region of the de 
Sitter universe with a periodic S 2 x R 1 spatial topology. 
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At this point we can discuss the computation of the classical energy term 

E SdS (M, b) = E dS (b) + AE§ S (M, b) ]elassieal 



(18) 



which can be computed by means of quasilocal energy. Quasilocal energy is defined as the value of the Hamiltonian 
that generates unit time translations orthogonal to the two-dimensional boundary, 



AElf (M, b) ]classical = ^ J d 2 x^ (k - k°) , 



(19) 



where | N\ = 1 at S 2 and k is the trace of the extrinsic curvature corresponding to the reference space, which in this 
case is the dS space. For practical purposes, however, it is convenient to embed both spaces (Sds and dS) into flat 
space and perform the subtraction procedure. To this purpose the radial coordinate x continuous on Ai is defined by 



dx = ± 



dr 



1 



2M(i 



(20) 



where the plus sign is relative to E+, while the minus sign is related to S_. The surfaces located at r + and ?* ++ 
are bifurcation surfaces denoted 5° and respectively. When M — 0, we obtain the embedding of dS space 

into flat space. In E + the evaluation of A£"£| s (M, b)i clas sical can be obtained as follows: first we consider the static 
Einstein-Rosen bridge associated to the SdS space fl6| , |l7|| 

- 2 - -N 2 (r) dt 2 + g xx dx 2 + r 2 (x) dfl 2 , (21) 



ds z 

where N, g xx , and r are functions of x defined by Eq.(|20|). Second, we consider the boundary S'+, located at 
x (r) = x + (R), and its associated normal = {h xx )^ 6^. The expression of the trace 



k = - 



7e(^"").„- 



gives for the SdS space 



k sds = 



r, x 



= -2 



VfV) 



2 / 2MG r 2 
\SdS r\ r b 2 



(22) 



(23) 



r \sds 

Note that if we make the identification iV 2 = 1 — 2M r G — p- , the line element (|2l]) reduces to the SdS metric written 
in another form. The same applies to the dS metric by putting M = 0. Nevertheless for our purposes the form of 
N (r) can be left unspecified. Thus the computation of E + gives 



AE$f 



(M,b) 



1 



\classical 



8ttG 



dn 2 r 2 



s 2 



-2v?W , 2 V /(r) i M =° 



(24) 



\r=R 



R 
G 



2MG R 2 / R 2 
1 R 52" _A ' 1_ ~ 



(25) 



where we have set M = in k sds to obtain the dS energy contribution. When R 3> b, AE^g S (M, b)^ classical ~ 
—iMb/R. Thus for every finite value of the boundary exceeding the cosmological radius, the classical energy acquires 
an imaginary component which will not be here considered^. In contrast, if we consider the approximation R/b -C 1 
and 2MG/R -C 1 @, we obtain 



2 To deal with this case it is better to introduce the quasilocal mass, defined as 

AM|f (M, b) lclassical = ^J 2 d 2 xNVt (k - fc°) 



R»b \ R \ b 
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AE§ S (M,b) lclaS8ical ^~ 



MG R 2 \ ( R 2 ^ 



1 R 2b 2 ) V 2b 2 



= M. (26) 



Thus the energy contribution limited to E + gives 

E SdS (M, b) = E dS (6) + AE s d dS (M, b)+ dassical = E dS (6) + M, (27) 

that is the dS space cannot decay into the SdS space because the associated boundary energy is different. This is 
in complete analogy with the Schwarzschild and the S-AdS cases. However if we look at the whole hypersurface 
E = S + U £_ the total classical energy becomes 

E sdS (M, b) = E dS (b) + E tot (M, b) 
= E dS (b) + AE s d dS (M, b)+ classical + AE s d dS (M, b)~ lassical (28) 

with 

AE S J S (M, b)+ las3ical = d 2 x^o (k - k°) , 

AE s d dS (M, b); classical =--^J s2 d 2 x^ (k - k°) . (29) 

Here the boundaries Si and S 2 . are located in the two disconnected regions Ai+ and M._ respectively with coordinate 
values x = x^ 1 and the trace of the extrinsic curvature in both regions is 



kSdS = , -2r, x jr onZ + 

1 2r, x jr on L_ v ' 



Thus one gets 



AE s d dS (M, b)f classical = { M M ° n m %_ , (31) 



where for E_ we have used the conventions relative to £_ and S_. Therefore for every value of the boundary _R, 
(provided we take symmetric boundary conditions with respect to the bifurcation surface), we have 

E SdS (M, b) = E dS (b) +M+ (-M) = E dS (b) , (32) 

namely the energy is conserved. As stressed in Ref. ]l6| , since we have a spacetime with a bifurcation surface, the 
quantities AE d g S (M, b)^ classical and AE d g S (M, &)|^; assica ; have the same relative sign, while the total energy is given 
by the sum AE d g S (M, b)^ classical + AE d g S (M, b)^ classica ^ The energies associated to the boundaries are symmetric 
and they have the same relative sign while the total energy reflects the orientation reversal of the two boundaries. 
Since the total classical energy is conserved we can discuss the existence of an instability. To this aim we refer to the 
variational approach to compute the energy density to one-loop [|l4|Jl5|Jl9| -^l| . 



III. ENERGY DENSITY CALCULATION IN SCHRODINGER REPRESENTATION 

In previous section we have fixed our attention to the classical part of Eq. (^|) . In this section, we apply the same 
calculation scheme of Refs. Jl4|,|l5| to compute one loop corrections to the classical SdS term. Like the Schwarzschild 
and the S-AdS case, there appear two classical constraints 



3 In Ref. |16| we have a subtraction instead of a sum. This is due to conventions adopted. 



5 



H 
W 



and two quantum constraints 



w = o 

W$ = 



(33) 



(34) 



for the Hamiltonian respectively, which are satisfied both by the SdS and dS metric on shell. H^ 1 = is known as 
the Wheeler-De Witt equation (WDW) . Our purpose is the computation of 



SdS 
£ 



H dS\ 
£ I 



(35) 



where H^. ds and i?^ 5 are the total Hamiltonians referred to the SdS and dS spacetimes respectively for the volume 
term |l4[j and ^ is a wave functional obtained following the usual WKB expansion of the WDW solution^. In this 
context, the approximated wave functional will be substituted by a trial wave functional of the gaussian form according 
to the variational approach we shall use to evaluate AE^g (M, b)^_ loop . To compute such a quantity we will consider 
on £ deviations from the SdS metric spatial section of the form, 



9ij = 9ij + h 



where 



9ij 



dx l dx 3 



, 2MG r- , , 

1 -pr ) dr 2 + r z dn 2 

r b z 



(36) 



(37) 



is the spatial SdS background. By setting M = in Eq.(g) on the same slice we will obtain perturbations also for the 
de Sitter metric. Thus the expansion of the three-scalar curvature / d 3 Xy/gR^ up to o (/i 2 ) gives 



/ d 3 xy/g 



-hAh + -h u Ah u - -^■■\ ! \,h'\ + ^hViVih li - ^Riah? + ^ ///.',, //'•' 



J d 3 xVd [\h H (RM - 6/6 2 ) h H - - A h u h u + \h (r^) h 



(38) 



where R^ is the three-scalar curvature on-shell. To explicitly make calculations, we need an orthogonal decomposition 
for both and hij to disentangle gauge modes from physical deformations. We define the inner product 

(h, k) := / ^/gG^hij (x) k kl (x) d 3 x 7 (39) 

by means of the inverse WDW metric Gijki , to have a metric on the space of deformations, i.e. a quadratic form on 
the tangent space at h, with 



Qijkl _ yikgjl + g ti g jk _ 2g i ig kl ). 
The inverse metric is defined on co-tangent space and it assumes the form 



(p. q) 



y/gGijkiP^ (x) q kl (x) d d x, 



so that 



G ijnm G nmkl = \ fas} + SfSi) 



(40) 



(41) 



(42) 



See also [E2H24] for other applications of the WKB approximation concerning black hole physics. 



G 



Note that in this scheme the "inverse metric" is actually the WDW metric defined on phase space. Now, we have the 
desired decomposition on the tangent space of 3-metric deformations |2qj2q]: 



1 



where the operator L maps & into symmetric tracefree tensors 



(/•CI,, = + V 3 -6 - -fly (V • . 



Then the inner product between three-geometries becomes 



(h, h) := / y/g&Mhij (as) ^ (x) d 3 z = 



(43) 



(44) 



With the orthogonal decomposition in hand we can define the trial wave functional 

* [ha (!>)] =AAex P (~ [{hK~ l h)^ y + {{L£)K- X {L0)[ y + (hK^h)^ 



Trace 

y 



(45) 



(46) 



where Af is a normalization factor. Since we are only interested in the perturbations of the physical degrees of freedom, 
we will only fix our attention on the TT (traceless and transverseless) tensor sector, therefore reducing the previous 
form into 



* [htj (H>)] = AAexp j-ij (hK- 1 h)^ v 



(47) 



This restriction is motivated by the fact that if an instability appears this will be in the physical sector referred to 
TT tensors, namely a nonconformal instability. This choice seems to be corroborated by the action decomposition of 
p7| , where only TT tensors contribute to the partition function^} To calculate the energy density, we need to know 
the action of some basic operators on VP [hij] ]l9|). The action of the operator h,j on \^) = \& [hij] is realized by 



hij{x)\y) = hij{H>)^[hij] 



while the action of the operator 7Ty on | , in general, is 

jry (a;) I*) = -i 



Shij (1?) 



*[hi 



The inner product is defined by the functional integration 

(*i I * 2 > = / [-Dhij]n{hij}^2{h k i} 



(48) 



(49) 



(50) 



and by applying previous functional integration rules, we obtain the expression of the one-loop-like Hamiltonian form 
for TT (traceless and transverseless) deformations [Il4 15 



H 1 



1 

Jm 



d 3 Xy /gG 



ijkl 



K 



-i± 



(x,x) ijkl + (A 2 ) a i K x (x,x) 



a f-L 
3 



The propagator K (x,x) iakl comes from a functional integration and it can be represented as 



ia kl 



h ia (^) Hi m 



N 



2X N (p) 



(51) 



(52) 



where h x a ("af) are the eigenfunctions of AjL and \n (p) are infinite variational parameters. 



3 See also [E8| for another point of view. 
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IV. THE SCHWARZSCHILD-DE SITTER METRIC SPIN 2 OPERATOR AND THE EVALUATION OF 

THE ENERGY DENSITY 



The Spin-two operator for the SdS metric is defined by 

(A 2 ); := -A^ Q + 2R a j - &/b 2 5 a } (53) 

where A is the curved Laplacian (Laplace-Beltrami operator) on a SdS background and is the mixed Ricci tensor 
whose components are: 

f IMG , , MG , 9 MG , ,1 
Rt = \ ^ + 2/& 2 ,-^+2/6 2 ,^ + 2/6 2 . (54) 



As stressed in the introduction, the form of the mixed Ricci tensor for the SdS space has the same dependence for 
the radial coordinate for both the Schwarzschild and the S-AdS spaces. Thus to evaluate the energy density, we are 
led to study the following eigenvalue equation 

(-A^ a + 2R a j - 6/b 2 5*) h\ = E 2 h) (55) 

where E 2 is the eigenvalue of the corresponding equation. In doing so, we follow Regge and Wheeler in analyzing the 
equation as modes of definite frequency, angular momentum and parity p9[ . The quantum number corresponding to 
the projection of the angular momentum on the z-axis will be set to zero. This choice will not alter the contribution to 
the total energy since we are dealing with a spherical symmetric problem. In this case, Regge- Wheeler decomposition 
shows that the even-parity three-dimensional perturbation is 



hr™ {r,$,<t>)=diag 



//,,., j i , r 2 K(r),r 2 sin 2 #K (r) 

r tr / 



Y m (■&,</>). (56) 



In this representation H (r) and K (r) behave as they were scalar fields. For a generic value of the angular momentum 
L, one gets 

' » G -£)H(r)=E 2 H H(r) 

(57) 

- 2 )K(r) = El K K(r) 

where Ef H and E 2 K are the eigenvalues for the H (r) field and the K (r) field respectively. The Laplacian restricted 
to E is 

A - ( i 2MG r2 } d2 + f 2r - 3MG 3 r \ d l Q> + 1 ) ( 58 ) 

\ r b 2 I dr 2 \ r 2 b 2 I dr r 2 




Defining reduced fields 



H(r) = ^; K{r) = k -^, (59) 



and passing to the proper geodesic distance from the throat of the bridge defined by Eq. ( po| ) , the system ( |57|) becomes^] 

-£jh (x) + (Vf (x) -$)h (x) = Efh (x) 

(60) 

-£k (x) + (V+ (x) - $) k (x) = Efh (x) 



with 



—j— ^ 1(1+1) 3MG , x 



The system is invariant in form if we make the minus choice in Eq. Q20 
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When r — ► ro > 



1(1 + 1) 3MG 



x (r) ~ ^2k+ (r - r+) V^ T (i) ► V 2 ; T — — = const, (62) 

where 

r 1 I If M ( r + ~ r -) ( r ++ ~ ? '+) f«o\ 

k+ = hm - \g 00 (r) = — 63) 

r^r + I Zb z r+ 

is the "inner" surface gravity associated with the smallest root. The solution of (|6(]) when r — ► ro > r + for both 
backgrounds is 

[2 

h (px) — k (px) — d — sin (px) . (64) 

This choice is dictated by the requirement that 

h (x) , k (x) — > when x (r) — > x (r+) ~ 0. (65) 

Thus the propagator becomes 

~L/ v V f°° 2 sm(px) sin(ra) Y l0 (#,<t>)Y VQ (■&,$) 

iv ± \ x iV) — 7T^2 / a PP — t\ t\ \ — \ 00 

2n 2 J r (x) r (y) A± (p) 

A± (p) is referred to the potential function Vj^ (x). Substituting Eq.(|66|) in Eq.([5l"|) one gets (after normalization in 
spin space and after a rescaling of the fields in such a way as to absorb l~>) 



E(M,b,\) = - s Y,Y, d pp 2 

1=0 i=l J0 



1=0 

where 



Ef(p,M,b,l) 

^ (p) + 



nJ , „ , 2 Z(Z + 1) 3MG 3 

El 2 (p,M,b,l)=p 2 + v 2 ; T ^^-_ (68) 

Ai (p) are variational parameters corresponding to the eigenvalues for a (graviton) spin-two particle in an external field 

and V is the volume of the system. By minimizing (p?]) with respect to Aj (p) one obtains A, (p) = [Ef (p, M, b, Z)] 2 
and 

E (M, b,X) = tt-jEE / *p a 2 A /£? (p, M, b, I) (69) 
;=0 i=l Ja 

with 



Z (Z + 1) 3MG 3 



For the SdS background we get 
where 



2 ' ■ V + c 2 _ + Jp 2 + 4 ) (70) 



2 _ Z (I + 1) 3MG 3 

C T ~ 12 T VS. ' 

'0 '0 

while when we refer to the dS space we put M = and c 2 = I (I + 1) /r 2 — 3/6 2 . Here the meaning of the value ro is 
that of maintaining the same boundary conditions to correctly compute the Casimir-likc energy. Then 
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y 00 r°° , 
77 1=0 •'° 



(71) 



Now, we are in position to compute the difference between (|70|) and (|71J). Since we are interested in the ultraviolet 
limit, we have 

AE (M, b) = E (M, b)-E (b) 



yV + c 2 _ + ^Jp 2 + c\ - 2yjp 2 + c 2 



—y 

4tt 2 ^ 



-n JO 



1=0 



dpp 3 



=Z) -2W1 
P 



(72) 



and for p >> c 2 ,c, we obtain 



V 

4^ 



00 poo 

£ / d pp 3 

1=0 Jo 



1 C 



1 c 



1 c 



I (2+ 

8 \p 



[p) + 4 \p 



2tt 2 8 



dp 

P ' 



where cf^ = 3MG/tq. We will use a cut-off A to keep under control the UV divergence 



dp 
P 



t dx , / A 
— ~ In 

11 a; \cm 



where A < m p . Thus AE (M, b) for high momenta becomes 



2tt 2 16 U^y 32tt 2 V r z J \3MG J 



and Eq.(J6l) to one loop is 



(73) 



(74) 



(75) 



{E SdS (M,b)-E dS (&)) 



r„ 3 A 2 



; r >r + 327r 2 ^ r 3 J 11 



X' 



3(r+r ++ ) (r+ + r ++ ) 



In 



2 (r 2 . + r+r++ + ri + ) rgA 



32^ 2 \2(r 2 + + r+r++ + r 2 + ) r% J \ 3 (r+r++) (r+ + r ++ ) 
where we have used Eq.(o3). On the other hand, when r — > r ++ 



xW.v^H Xf (x) — ^±11 T 3MG 



const, 



++ 



++ 



where 



(76) 



(77) 



(78) 



is the "outer" surface gravity associated with the largest root. By repeating the steps going from Eq.(|67|) to Eq.(|7q) 
with 
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E 2 2 (p,M,b,l)=p 2 + \ T — 3 s , (79) 



we obtain 



(*» (M, », - (6,)^ + = (^) (^) 



V / 3r+(r + +r++) V, /2 (r| + r+r++ + r\ + ) r^A 2 ' 



32tt 2 ^ 2 (r 2 + r+r ++ + r\ + ) r\ + J \ 3r+ (r+ + r++) 

Like the Schwarzschild and the S-AdS cases, we observe that 

lim lim AE (M, b) ^ lim lim AE (M, b) . (81) 

M— >0r — >r + ro— *T+ M— »0 

This behavior seems to confirm that quantum effects come into play when we try to reach the inner throat. Differently 
to the Schwarzschild and S-AdS cases, here we have two scales corresponding to the two throats (horizons) of the 
metric. This is an artifact of the approximation we have adopted to deal with the differential equations of the system 
(|60|). By defining a scale variable x = 3MG/ (r^A 2 ) and a scale variable y — 3MG/ (r^_ + A 2 ), we obtain 

E SdS (M, b) - E dS (6) = {E SdS (M, b) ~ E dS (b))^ ro>r+ + {E SdS (M, b) - E ds (6)) r=ir++ 



V 
"32^ 



/ 3MG\ fr 3 Q A 2 \ {3MG\\ /r^ + A 2 



VA 4 

= —^[x 2 \nx + y 2 \ny\=AE{x,y). (82) 

A stationary point is reached for x = y = 0, namely the dS space and another stationary point is in x = y = e~2. 
This last one represents a minimum of AE{x 1 y). This means that there is a probability that the dS spacetime will 
be subjected to a topology change and it will produce a SdS wormhole with a black hole pairQ generated on the 
hypersurface S. To see if this is really possible, we have to establish if there exist unstable modes. 

V. SEARCHING FOR NEGATIVE MODES 

In this paragraph we look for negative modes of the eigenvalue equation (|53]). To this purpose we restrict the 
analysis to the S wave. Indeed, in this state the centrifugal term is absent and this gives the function V (x) a potential 
well form, which is different when the angular momentum I > 1. Moreover the potential well appears only for the H 
component, whose eigenvalue equation is 

-A - — — - ^ J H (r) = -E 2 H (r) . (83) 
A is the operator A; of Eq.(|58|) with I = and Eq > 0. By defining the reduced field h(r) = H (r) r, Eq.(|83j) becomes 



, i _2MG_ph\(-3MG + - .. 



dr \ V r b 2 dr 



^ _ 2MG 



7 See also [Eol for an approach of black hole pair creation based on the effective action. 
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where E 2 — —3/b 2 + E 2 . By means of Eq.(§0), one gets 



dx d 
dr dx 



2MG r 2 dh dx 
r b 2 dx dr 



3MG - 2 



1 - 



2MG r> 



Near the throat 



where we have used Eq. 
obtain p ~ 1 + y 2 where 



Eq.([S5|) becomes 



d Uh\ 
dx \dx J 



3MG 



E 2 \h = 0. 



x (r) ~ ± 



— -i , 



(85) 



(86) 



/«+ v v r + 

defining the inner surface gravity. By defining the dimensionless variable p = r/r+, we 



y = v / K+a;/y / 2r + = k + x. 



(87) 



d f dh 



dy \dy 



3MG 

r+P 3 (y) 



E 2 h 



d 2 h 
dy 2 



3MG 



k\r\ (1 + y 2 f 



+ A \h = 0, 



where A = E jk Expanding the potential around y = 0, one gets 



cPh 
dy : 

~dy^ 



+ 



-^(l-3y 2 ) + X)h 



2 2 3MG , . , 



(89) 



(90) 



with w = J9MG/ (k+t+). In this approximation we have obtained the equation of a quantum harmonic oscillator 
equation whose spectrum is E n = fuo (n+ 4). Since we are using natural units, h is set to one and 



A„ = 3MG/ (R 2 + r\) - yj9MG/ (k 2 + r 3 + ) (n+^ 



After some algebraic calculation, with the help of relation fll4|), we obtain 
An = 3 



2r++ (r+ + r+H 



2r ++ (r+ + r + 



1 I 1 

tiH — 

(2r+ + r++) (r ++ - r+) IV (2r+ + r ++ ) (r ++ - r+) ^2 V 2 



Let us examine the first two eigenvalues. The first one is 
A = 3i 



2r++ (r+ +r ++ ) 



2?-++ (r+ +r ++ ) 

(2r+ + r ++ ) (r++ - r+) \]j (2r+ + r++) (r ++ - r+) 

Since the eigenvalue must be positive, the following inequality must hold 

1 



' 2r ++ (r++r ++ ) 1 
(2r + + r ++ )(r ++ -r+) 2 



7 

4 + 4 



7 2 

4 + " 



<f + > °< 



(91) 



(92) 



(93) 



(94) 
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which is verified V6> S [§ ,tt). To proof that there is only one eigenvalue, we look at the second one 




2r++ (r+ +r++) 
(2r + +r ++ ) (r ++ - r+) 



(95) 



This implies the inequality 



18r£ > r+r ++ + r\ + , (96) 
which is verified when 9 £ (2.1708, 7r). Thus we have only one eigenvalue when 6* G [l^-, 2.1708). In terms of E 2 we get 

^ = _ 3/6 2 _ 3r++ ^++ r++) + JLf V2r ++ (r + + r ++ ) (2r + + r ++ ) (r ++ - r + ). (97) 

Note that if we repeat the same calculation for the outer throat, i.e. we discover the absence of negative 

eigenvalues. Thus we can conclude that there is only one eigenvalue with the restriction that 9 € [^,2.1708). 
According to Coleman j^], this is a signal of a transition from a false vacuum to a true one. It is evident that when 
we consider the limit where r ++ — > r+, an infinite number of eigenvalues enter in the discrete spectrum. This is a 
consequence of the approximated parabolic potential which is used to describe this extreme situation. To better deal 
with this problem, we introduce the Nariai metric. 



VI. THE NARIAI METRIC SPIN 2 OPERATOR AND THE EVALUATION OF THE ENERGY DENSITY 

When r + = r ++ , the metric becomes degenerate and the function / (r) of Eq. (|llj) becomes 

fe(r) = -^(r-f) 2 (r + 2f), (98) 

where f = b/y/3. Since f e (r) < everywhere r becomes a time coordinate and t becomes spatial. Nevertheless, this 
is an artifact of a poor coordinate choice. To see what happens, we follow Ref. Q and we let 9M 2 G 2 A = 1 — 3e 2 so 
that the limit r + — > r ++ corresponds to e — > 0. We define a new radial coordinate 6\, and a new time coordinate 4>\, 
by 

COS0! = ^( r -f) fc = t^, (99) 

be be 

such that r + = f — e\f%jb and r ++ — f + eV3/b. To first order in e the metric assumes the form 

b 2 ( 2 \ b 2 ( 2 \ 

ds 2 = -— ( 1 + -ecos6»i j sin 2 6\d<f)\ + — I 1 - -£cos6»i j d9\ 

b 2 

+ — (1 - 2ecos9 1 )dn 2 , (100) 
3 

describing a nearly degenerate SdS metric with two distinct roots. The related surface gravities, to first order in e 
assume the expressions || 

«+,++ = 4= ( 1= F^e cos O^j , (101) 

where the upper (lower) sign is for the cosmological (black hole) horizon. When e — > 0, the line element is 

b 2 

ds 2 = — (- sin 2 0id</> 2 + d6» 2 + ci^ 2 ) (102) 
and the surface gravities degenerate into one with the value 



K = 7T (103) 
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This is the Nariai metric Q with topology H 2 x S 2 . Its Euclidean form has the well known topology S 2 x S 2 and 
form 

b 2 

d s 2 = — (dtt\ + dfl 2 ) . (104) 
3 

Every constant time section has the topology S 1 x S 2 with the throats of the same size as illustrated in Fig.|| 



A 



r=i>r ++ r=r=r„ 

FIG. 3. The geometry of the constant time slice degenerate Schwarzschild-de Sitter (Nariai) spacetime with a polar angle 
suppressed. Isometric copies of this surface can be smoothly joined at the throats, producing a periodic S 2 x S 1 spatial topology. 

However we are interested for the Lorentzian version. In this case, Regge- Wheeler decomposition shows that the 
even-parity three-dimensional perturbation is 

h% m (0i, 0, 4>) = jdiag [H (0x) , K (6 1 ) , sin 2 OK (6 1 )] Y l0 (0, 0) . (105) 

For a generic value of the angular momentum L, one gets 

£(~d 2 i +l(l + l)-2)H(e 1 ) = E 2 H(6 1 ) 

(106) 

£(~d 2 i +l(l + l))K(9 1 ) = E 2 K(9 1 ) 
where the associated three dimensional mixed Ricci tensor is 

R i = {°44}- (107) 

Also in this case the unstable mode appears for the I = case leading to the eigenvalue equation 

(-9 2 1 -2) J ff(0 1 ) = -Ai/(0 1 ), (108) 
where A = b 2 E 2 /3. The eigenvalue is easily determined and its value is 

A = 2, (109) 

with eigcnfunction 

H (0i) — const. (110) 

Since the range of integration is finite due to the periodicity of the argument, the eigenfunction is normalizable. For 
the other component we have no solutions at all, because the operator is the same of a free particle having only a 
continuous spectrum. For completeness, we calculate the energy difference between the SdS background in the Nariai 
form and the dS metric in the stable sector. The total energy in the presence of the Nariai metric is 



, , + «,„ , ,L, . 31(1 + 1) 



E Nariai (b) = ^-yj dpp 2 \\lp 2 + ^ yr ^-6/b 2 + ] Jp 2 + ^ yr ^ L (111) 



while for the pure de Sitter metric, we have 



VI f 

;_n J U 



x '- x / ; 3/ (/ + :i ) 



The Casimir-like energy becomes 



p 2 + y b2 ' - 3/b 2 . (112) 
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AE (b) = E Nariai (b) - E dS (b) 



V 1^ f°° , a 

1=0 J ° 



9 31(1 + 1) e/ ,„ / , 31(1 + 1) 
j 1 i - 6/6 2 + Wp 2 + — 



b 2 



b 2 



-2\/P 



, 31(1 + 1) 
6 2 



3/6 2 



Whenp 2 >> 6/6 2 we write 



Introducing an UV cut-off one gets 



and AE (b) for high momenta becomes 



AE (6) 



V 9 f°° dp 



2tt 2 4& 4 7 p 



dp f A b dx 
o P Jo x 



In (A 2 b 2 



AE(b) --J-\n(A 2 b 2 ) = 

V ' 2ir 2 16b 4 V ' 



V 9 
32tt 2 fr* 



In (A 2 6 2 



(113) 



(114) 



(115) 



(116) 



VII. BOUNDARY REDUCTION AND STABILITY 



An equivalent approach to Eq.(p4[) can be set up by means of a variational procedure applied on a functional whose 
minimum represents the solution of the problem. Let us define 



X 

j(h,E 2 ) =1 / 'dx 



(dh{x)\ 3MG 
\ dx 



3(x) h2{x) 



E 2 

+ — I h 2 (x) dx, 



where dx is given by Eq.(||). Eq.(|8J) is equivalent to find the minimum of 



E 2 = °- 



fdx 



/ h 2 (x) dx 
o 



For future purposes, we use the boundary conditions 

h(x) =0 

When 2MG/r < 1 Eq.@ becomes 

d 2 h 



dx 2 



E 2 h = 



and the solution is 



h (x) = Aexp y—Exj + Bexp LExj 



(117) 



(118) 



(119) 



(120) 



In the approximation of Eq. (|ll!j) , we can suppose x so large that the increasing exponential has to be eliminated in 
such a way to take only h (x) — Aexp (—Ex \ . If we change the variables in a dimensionless form like Eq.(|S7j), we get 
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M = — 

K 



y 

fdy 

o 



(ihM) 2 - JMG (y) 
\ dy J r^n^p 3 (y) W 



y( a ) 

J dyh 2 (y) 



(121) 



The asymptotic behaviour of h (x) suggests to choose h (A, y) = exp (—Ay) as a trial function, and Eq.(121) becomes 



M (A) = A 



Close to the throat exp (—2Xy) ~ 1 — 2Ay and 

M (A) = A 2 



f 4h exp (-2 Ay) 



3 £.2 l-cxp(-2A;/) 



2A 



UIG 9MG 



V , -2 
2A +2/ . 



(122) 



(123) 



The minimum of p (A) is reached for A = f jpr|i y^j 3 assuming therefore the value 



where 



p(X)=3[ J Dy) -3D + 3Dy 2 , 
2r ++ ( r ++ r ++) 



(124) 



D 



MG 



k\ (2r+ + r ++ ) (r++ - r+) 



2oob(D (cos(^)+cos(f)) 



(2cos(^)+cos(f))(cos(f) 



cos I 



")) 



(125) 



If we consider the value of 9 such that only one eigenvalue appears, we obtain D = 2.25 and Eq.(124)becomes zero 
for y c — .463 68 corresponding to p c = 1.215. This means that the unstable mode persists until the boundary radius 
p falls below p c . 



VIII. SUMMARY AND CONCLUSIONS 

In this paper we have extended the computation of the Casimir-like energy to the case of the Schwarzschild-de Sitter 
(SdS) background with the de Sitter (dS) space as a reference space. This evaluation has been done to one-loop in 
the TT (transverse, traceless) sector which is the gauge invariant part of the quantum fluctuation of the gravitational 
field. As stressed in the introduction to correctly compute the Casimir energy we need to subtract field configurations 
which have the same asymptotic properties and the same asymptotic boundary conditions. For the SdS metric and 
the dS metric, this is the case. In this context a lot of work has been done; for example in Refs. [0-0], it has been 
shown the existence of one negative mode in the TT sector when the saddle point approximation is considered: a 
clear sign of instability. In particular, the instability has been related to the probability of creating a black hole pair 
|^,|lj. However this particular result has been obtained by looking at the partition function and therefore with the 
introduction of an equilibrium temperature, that in the case of the pure Schwarzschild and flat metric can be imposed 
to be equal, but in the present case (i.e. the SdS metric and the dS metric) it cannot. Therefore as stressed in Ref. 
||, it is not very clear how these spaces having a different temperature (periodicity) can be compared. On the other 
hand if we adopt the hamiltonian approach we can avoid the introduction of a temperature and it is possible to build 
a scheme where the classical contribution is conserved; this point is fundamental to discuss the instability g]. What 
we have found, in our framework, is the well known existence of an unstable mode in the S wave for the extreme 
SdS metric (Nariai)R but we have also found that an unstable mode exists also for the SdS metric. To interpret such 



s For this last point a discussion can be found in Ref. 
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an instability as a decay process nucleating a black hole pair, we need to show the existence of only "one negative 
mode" @]. Unfortunately, in our approximation we have discovered a dependence of the number of negative modes 
on 9 the variable defined in Eq.(|l2|). In particular, we have seen that the closest is the approach to the extreme 
value 9 — it, the highest is the number of negative eigenvalues falling into the negative spectrum. This abundance 
of negative eigenvalues can be limited if we restrict the variability of 9 into the interval [5,2.1708); in this range 
Coleman arguments can be applied. However the proliferation of negative modes in the nearly extremal SdS metric is a 
consequence of th e ap proximated potential of the eigenvalue equation. Indeed by using nearly degenerate coordinates 
like those in Eq.( |l00| ), only one negative eigenvalue appears. It is interesting to observe that the appearance of a 
negative mode, even for the SdS metric, can be related with the production of a sub-maximal black hole pair [ p2[ . 
What is interesting to observe is the existence of a critical radius p c below which the instability disappears. This 
could open the possibility of an existing foam-like space composed by copies of bubbles in analogy with the model 
discussed in Ref. |53|. Even in this case, a dependence of the ultraviolet cut-off is present. This is principally due to 
the non renormalizability of quantum gravity. However, the fact that the same divergent behaviour appears also in 
this case it is a signal of a more general situation typically concerning the spherically symmetric metrics, that is every 
spherically symmetric metric describing a wormhole (black hole) can produce in the spectrum of quantum fluctuations 
in a semiclassical approximation one negative mode, provided the boundary conditions be energy conserving. 
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APPENDIX A: KRUSKAL-SZEKERES COORDINATES FOR SDS SPACETIME 



We have defined the SdS line element in Eq.(^). To introduce the Kruskal-Szekeres [34 -36 type coordinates we 
consider the following transformation 



ds 2 



2MG 



r 

&2 



[dt 2 - dr* 2 ] + r 2 dQ 2 



2MG 



— I dvdu + r 2 (u, v) dfl 2 
r b z 



(Al) 



where v = t + r* is the ingoing radial null coordinate and u — t — r* is the outgoing radial null coordinate. The 
"tortoise coordinate" r* is defined by 



dr* = -- 



rb 2 dr 



(r - r+) (r - r ++ ) (r + r+ + r ++ ) 
and by means of the surface gravity associated to each root, we write 

11 1 1 



/ 2n + (r-r+) 2n ++ (r - r ++ ) 2k_ (r + r+ + r++) ' 



(A2) 



(A3) 



where / (r) has been defined by Eq. (g) , while k ++ and k + by Eq. ( |78D and Eq. (| 
gravity associated with the negative root is 

(2r+ + r ++ ) (2r ++ + r+) 



2b 2 (r+ +r ++ ) 



respectively. The last surface 

(A4) 



Thus 





r 










1 




r ^ 




2k+ 


r+ 




r++ 





1 



(A5) 



To avoid singularities we can define Kruskal-Szekeres type coordinates 



V 



++ _ 



exp k++v 



U = — exp— k ++ u. 



(A6) 
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These coordinates do not cover r < r + because of the coordinate singularity at r = r + (and U ++ V ++ is complex for 
r < r_|_), but r — r + and a similar four regions are covered by the (U + , V + ) Kruskal-Szekeres-type coordinates to this 
case. 

V + — expn + v U + — — exp— n + u. (A7) 

For the H — h sign we have 

+ / r A K ~ l r 



U++V ++ = -exp (v-u)) = -exp (2 K++ r*) =- ( — - lj I - — + 1] j— - t) (A8) 
and the respective line element is 



^,=- (r++ r+ +)r+r++ ( + lV " f — - 1 V " + df/+W+++r 2 ^++.r- 'j./o^ 



i-2±±. / \ 1 

r 



1 — -1 d(7++^+ + + r 2 (C/++,^++)dr! 2 , (A9) 



2MG 

while for the + sign we have 

= - exp (k + (u - u)) = - exp (2/c+r*) = - (r - r + ) (r + r + + r ++ ) ~- (r ++ - r) K ++ (A10) 
and the associated line element is 

= _ (r + +r ++ )r + r ++ / \ 1 / _ _r\ ^ ^ ^ + ^ 



= f + lV K ~ (l - —) 1 K++ dU+dV+ + r 2 (U+,V+) dQ 2 . (All) 

K+r \r++r ++ J \ r ++ J 

The conformal Penrose diagram of the SdS space is shown in Fig.|]. Regions limited by r + < r < r ++ lie between 
the black hole and cosmological horizon. Regions r > r ++ correspond to an asymptotic de Sitter region and region 
r + > r to the black hole interior. 



r=0 




FIG. 4. Penrose diagram for the Schwarzschild-de Sitter spacetime. 
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